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Abstract 



In addition to the Riemannian metricization of the thermodynamic state 
space, local relaxation times offer a natural time scale, too. Generalizing 
existing proposals, we relate thermodynamic time scale to the standard ki- 
netic coefficients of irreversible thermodynamics. Criteria for minimum en- 
tropy production in slow, slightly irreversible processes are discussed. Euler- 
Lagrange equations are derived for optimum thermodynamic control for fixed 
clock-time period as well as for fixed thermodynamic time period. Only this 
latter requires constant thermodynamic speed as the optimum control pro- 
posed earlier. An easy-to-implement stepwise algorithm is constructed to re- 
alize control at constant thermodynamic speed. Since thermodynamic time 
is shown to correspond to the number of steps, thus the sophisticated task of 
determining thermodynamic time in real control problems can be substituted 
by measuring ordinary intensive variables. Most remarkably, optimum paths 
are Riemannian geodesies which would not be the case had we used ordinary 
time. 



1 Introduction 



Standard equations of irreversible thermodynamics have been known for 
many decades. Investigations on optimizing finite-time thermodynamic pro- 
cesses controlled by external reservoirs date from the 80 's. An instructive 
exposition of the problem at the time was published in Ref. Basically, 
the problem consists of finding the best path in the state space along which 
one drives the system from a given equilibrium state to another. 

A particular approach to the problem of optimally controlling finite-time 
thermodynamic processes takes its origin from the natural geometric struc- 
ture of the thermodynamic state space [0, |3|, f|. It would seem straightfor- 
ward to expect that the geodesic path should somehow be related to the 
optimum path connecting the given initial and final equilibrium states. A 
few years ago it was pointed out 0, |(J that the natural relaxation time r plays 
a fundamental role in devising optimum cooling strategies, e.g. in computer 
simulated annealing. Quite recently, Andresen and Gordon J7J have shown 
that the strategy of constant thermodynamic speed || is related to a certain 
minimum of entropy production. In the present work we reconsider these 
ideas and make a definite progress. 

The problem itself is illustrated by the finite-time cooling process 
(Sect. 2), its dynamics is described by a phenomenological cooling-equation. 
Then, applying kinetic equations from standard irreversible thermodynamics, 
we generalize the concept of thermodynamic time and speed for any num- 
ber of control variables (Sect. 3). Furthermore we invoke standard Euler- 
Lagrange equations to obtain the finite-time path of minimum entropy pro- 
duction and we make explicit the role of variational conditions. Most impor- 
tantly, we prove that the optimum path is geodesical so a broader generaliza- 
tion of the old principle of constant thermodynamic speed is achieved. Even 
an older belief about the significance of Riemannian metric in thermodynam- 
ics might get justified (Sect. 4). Finally, we construct an iterative algorithm 
to realize thermodynamic processes at constant thermodynamic speed, also 
giving a genuine control-theoretic interpretation of thermodynamic time itself 
(Sect. 5). 
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2 Example: Finite-time cooling 



Consider a system in thermal contact with a large reservoir and let T and T r 
denote their respective temperatures. In general, both T and T r will depend 
on time t. At the initial time ti, assume the system and the reservoir are in 
equilibrium with each other at temperature Ti, i.e., 

T r (ti) = T(ti) = Ti. (1) 

Then start to decrease the reservoir's temperature T r from Tj to Tf, con- 
suming fixed finite time tf — ti, i.e., choose a given function T r (t) so that 
T r (tf) = Tf. The system's temperature T(t) decreases due to permanent 
heat transfer to the reservoir and will always be retarded with respect to the 
reservoir's current temperature T r (t) by some positive AT(t) = T(t) — T r (t). 
Throughout this paper, we consider slightly irreversible processes when, e.g., 
T r (t) changes slowly enough to allow the heat transfer to satisfy Newton's 
law and the following equation is expected to drive the permanent relaxation 
of the system's temperature T: 

provided the local relaxation time r changes little between T and T r . It will 
be useful to re-scale the clock-time parameter. Introducing thermodynamic 
time £ was proposed earlier in Refs. 0, |6], [7J • The relation of the two scales 
relies upon the local relaxation time r(T(t)) along the cooling process: 

d£ = dt/r. (3) 

In the new variable, Eq. (|2|) takes a simple form: 

T'.f = -AT. (4) 

We note that this equation has the following explicit solution with the initial 
condition & = 0: 

T(0 = e-«(T i + ^T r (£')e r O- (5) 

Now the basic goal is to single out "optimum" cooling paths. Following 
e.g. Andresen and Gordon 0] one requires that the optimum cooling happen 
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with maximum reversibility, i.e., at minimum total entropy production. The 
entropy production rate of the cooling process is 

S = C (T)f(i-l) (6) 

where C{T) is the specific heat of the system. In case of sufficiently slow 
cooling this expression reduces to 

S = C(T)(f/T) 2 T(T) (7) 

where we applied Eq. (|2|). The common criterium of optimum is 

Sdt = min. (8) 

u 

where the overall time £/— i, of the process is fixed. This optimum is achieved 
when 

S = ^pT 2 r(T) = const. (9) 

There is, however, a remarkable alternative to this optimum, because one 
can choose different boundary conditions. Instead of clock-time £/ — ij the 
thermodynamic lapse £/■ — £j of the cooling can be fixed as well. Then the 
optimum cooling becomes different; it will correspond to constant entropy 
rate versus thermodynamic time: 

S' = tS = ^p-T' 2 = const. (10) 

This condition has a challenging geometrical interpretation |7J. 5" is the 
square of the thermodynamic speed of the cooling process: 

S'=\\T'\\ 2 , (11) 

calculated with the entropic metric defined by the quadratic norm 

WdTf = ^P(rfT) 2 . (12) 

Hence the corresponding optimum process is called cooling at constant ther- 
modynamic speed (cf. Refs. [|, || 0). In Sect. 4 we shall prove that the 
principle of constant thermodynamic speed also applies for optimum finite- 
time thermodynamic processes affecting more (than one) variables. First, in 
Sect. 3 we must generalize the concept of thermodynamic length, time, and 
speed for such finite-time processes. 
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3 Thermodynamic length and time 



Let the equilibrium states of a given thermodynamic system be character- 
ized by the n + 1 extensive variables; the vector X = (X 1 , X 2 , . . . , X n ) will 
parametrize the manifold of state space while X n+1 remains fixed. Following 
Refs. [fj, |3|, |J one defines a metric tensor g on the manifold of equilibrium 
states, derived from the entropy S(X): 

d 2 S(X) 

»w = -mU- (13) 

If one took the entropic intensive variables Yk = dS(X)/dX k ; k = 1,2, ... ,n 
instead of the extensive variables X then the metric tensor would be the 
inverse g^ 1 = [g lh ] of g = [g^]. Hence, in obvious notations, the thermody- 
namic line-element square takes the following alternative forms: 

\\dX\\ 2 = (dX\g\dX) = \\dY\\ 2 = (dYlg-^dY). (14) 

Consider now a certain path {X(t); U < t < tf} in the thermodynamic state 
space. The thermodynamic length I of the path takes the (alternative) forms 



dx\\ = / ||dy||. (15) 



Obviously, the length of a path is independent of the choice of coordinates and 
even of the parametrization of the path itself. No intrinsic relation manifests 
itself between the thermodynamic length and the clock-time parameter t. 

In order to obtain a natural time scale along a given path, we first have to 
invoke standard concepts of irreversible thermodynamics. Consider a reser- 
voir in equilibrium at some state variables X r and bring it into contact with 
the system which is in equilibrium at X. Then, the state of the system will 
converge to the state of the reservoir. The standard form of the relaxation 
equations reads |§: 

^=7^ AF fc (16) 

where AY = Y r — Y is the deviation from the equilibrium in terms of the 
entropic intensive variables, and 7 = [y lk ] is the matrix of kinetic coefficients; 
it is symmetric and positive. During this process of relaxation, entropy S 
will be produced at rate 

S = (X\AY). (17) 
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If AF is small, one can write 



AY = -gAX. (18) 

The relaxation equation ( |TBD takes the alternative forms: 

X = -jgAX or Y = -gjAY. (19) 

Using the equation (|1|) and the equation 

Y = -gX (20) 

alternative expressions of entropy production rate (|17D follow: 

S=(X\ 1 ~'\X) = (Y\(g 1 g)-'\Y). (21) 

Based on the above standard equations, an effective relaxation time r, 
depending on the path's local direction X, can be defined. Recall from Sect. 2 
that the system's path X(t) is driven by a certain reservoir path X r (t) and 
the system's retardation AX = X— X r is proportional to the current velocity 
X of the process: 

AX = -^gy'X (22) 

In fact, this equation is formally identical with the relaxation equation (|19|) 
which remains valid if X r becomes a slowly varying function of time. Observe 
that the longitudinal (i.e. parallel to X) projection of the above equation 
implies a certain effective relaxation time r. Introduce the longitudinal com- 
ponent of the retardation: 

AX" = -L(AXMX) (23) 
II 

Substituting Eq. (p2|) yields the longitudinal relaxation equation (note that 
Xll = ||X||): 

AX 11 = -r||X|| (24) 

where 

' (AX|<7|X). (25) 



IX 
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This can be rewritten in the following equivalent form: 



t=-- i -(X\j-i\X). (26) 

Indeed, as it is clearly seen from this form, the effective relaxation time r is 
positive and only depends on the direction of X but not on its magnitude 
(as long as it is moderate). Eq. (|2~6"D have the compact form 



r = (27) 

showing up r's invariance if one changes the representation of the states from 
extensive variables X to intensive ones Y, for instance. 

Having introduced the effective relaxation time r, the notion of thermody- 
namic time £ can now be extended for paths in more dimensions. Formally, 
we retain the old definition @ of the element of thermodynamic time £: 

d£ = dt/r (28) 

which now depends on the direction of the speed X. Sometimes, a vector of 
thermodynamic speed X' (or Y') will be more useful than X (or Y): 

. dX . dY 

X' = — = tX or Y' = — = tY. 29 

d£ d£ v ' 

The common (scalar) thermodynamic speed || corresponds to the invariant 
norm(s) ||X'|| = of the vector(s) (|29|) so an extension of the notion of 
thermodynamic speed for more dimensions has been performed. The entropy 
production S' = tS per unit thermodynamic time can also be considered. 
From Eq. ( p7|) we obtain: 



S' = tS= \\X'\\ 2 (30) 

This means that the dimensionless entropy production rate S' is equal to the 
squared thermodynamic speed, as is expected from the single variable case 
in Sect. 2. 
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4 Optimum paths minimizing entropy pro- 
duction 



Consider a certain path {Y(t)} corresponding to a finite-time thermodynamic 
process connecting the initial state Yi = Y(ti) with the final one Yf = Y(tf). 
Remember that a thermodynamic path {Y(t)} is the solution to the "cooling" 
equation (^|) having now the following general form [cf. Eq. (|19|)1: 

Y = 91 (Y r -Y), (31) 

driven by the given reservoir path {Y r (t)}. The system's path has a small 
retardation AY(t) behind the reservoir path. In the slow process approxi- 
mation Eq. ([2lD applies and the total entropy production will depend on the 
path as follows: 

Sfi = f tf Sdt = f 1 '(YKg^y^dt. (32) 

Let us find the path minimizing the overall entropy production, among paths 
connecting the fixed initial and ending points at a fixed time lapse £/ — U. 
An analogy with Lagrange's variational principle can be established if we 
identify the Lagrange-function as 5/2. Then, minimizing paths are found to 
obey the following Euler-Lagrange- equations: 

= (33, 

Obviously, the entropy production rate is an integral of the above dif- 
ferential equation: 

S = (YKg^g)- 1 ^) = const. (34) 

which is the mathematical counterpart of the energy conservation rule in 
mechanics. In our case this means that the entropy minimizing path corre- 
sponds to a constant entropy production rate, as was pointed out in Sect. 2 
for the single variable cooling process. 

An equivalent Euler equation could be obtained had we chosen the ex- 
tensive variables X to parametrize the paths. Actually we have chosen the 
intensive ones since in typical experimental situations the reservoir's intensive 
variables are under our control (cf. Sect. 2). 
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In certain cases (see Sect. 5) it would be interesting to find the path 
minimizing Sfi at the condition that the thermodynamic time — £j = 
f£f dt/r be kept fixed and the clock-time tf — ti might be varied. In this case, 
one suitably replaces all t-dependences by ^-dependences. From rate (|30| ). 
one obtains the following equation for the total entropy production: 



S fi = f I S'd£ = I^WY'fdi. (35) 



At fixed the minimum of entropy production is achieved if the path 

Y(£) satisfies the following Euler-Lagrange equations: 

(36) 

The entropy production rate versus thermodynamic time £ will be an integral 
of this Euler-Lagrange equation: 

S' = || F' || 2 = const. (37) 

Thus the minimizing path corresponds to constant S'. On the other hand, 
S' is equal to the squared invariant norm of the thermodynamic speed Y' 
defined by Eq. (p9|). Hence the optimum path is of constant thermodynamic 
speed. An equivalent result could be obtained in A-variables. The constancy 
of the thermodynamic speed is merely a consequence of a much remarkable 
feature of optimum paths: they are geodesies of the Riemann-metricized 
manifold of thermodynamic states. 



5 An easy control of optimum cooling 

It follows from the previous Section that the Riemannian geometric structure 
of the thermodynamic state space has an intrinsic relation with the processes 
of maximum reversibility at fixed thermodynamic lapse rather than clock- 
time. We are going to show that optimum cooling processes can easily be 
controlled by using thermometers instead of clocks! The method is straight- 
forward to implement for the simultaneous control of more thermodynamic 
parameters. 
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We propose a stepwise control strategy to approximate the optimum pro- 
cess at constant thermodynamic speed. Assume that we change the reservoir 
temperature stepwise between T and Tf\ 

Ti = To = T t .q > T r i > T r 2 > . . . > T t n = Tf, 

and the cooling strategy goes like this. First we lower the reservoir tem- 
perature to T r i suddenly, and we wait for the system to become equili- 
brated to a given extent e, i.e. its temperature T will be as close to T rl 
as Ti = eT rl + (1 — e)Tj. Then we move to the next iteration by lowering 
the reservoir temperature to T r 2 and letting the system's temperature T to 
equilibrate to the same extent e, identified by T 2 = eT r2 + (1 — e)T\. In 
general: 

T n+1 = eT r , n+l + (1 - e)T n , (n = 0, 1, . . . , N - 1). (38) 

If all steps are so small that the change of the relaxation time t(T) is negli- 
gible at the step's scale then the relaxation equation (H) holds and each step 
will have the same lapse 

*( = log(l-r' (39) 
of thermodynamic time. Hence, the number of small steps required to 
realize the cooling process from Tj to Tf, provided the quality of each relax- 
ation has had the common characteristic value e, will be proportional to the 
thermodynamic time lapse £f — & of the cooling path: 

N = ^ - - . (40) 
|log(l-e)| 1 ] 

The smaller the defect of relaxations 1 — e the bigger number of iterations 
will be necessary to achieve the same cooling Tj — > Tf. 

To assure constant thermodynamic speed in average we choose steps 
T rl , T r2 , ... in such a way that the same thermodynamic length 5£ defined by 



(|i~2D belong to all corresponding segments ||Tj — Ti||, ||Ti — T 2 ||, . . .: 

T n+1 - T n | 2 = 5f, (n = 0, 1, . . . , N - 1) (41) 



J>2 



From equations (p8| , ^l|) one obtains 

We 
/CXT 



T r , n+1 = |1- S J^— \T n , (n = 0,l,...,Ar-l). (42) 
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Now we can summarize the iterative algorithm of cooling. Given the cur- 
rent temperature T n of the system, one decreases the reservoir's temperature 
according to Eq. (|3J) and lets the system relax until the condition fl3g|) be- 
comes valid. This procedure is then iterated for T n+ %, e.t.c. One can see by 
inspection that the same average thermodynamic speed 

" = \^h\ (43) 

belongs to each steps and consequently to the whole process as well. So this 
stepwise algorithm approaches the theoretical optimum process of constant 
thermodynamic speed ||. We need no clocks to measure time but thermome- 
ters to measure temperatures and the a priori knowledge of the specific heat 
function C(T) of the system in hand. 

In the general case when more thermodynamic variables are to be con- 
trolled one has to derive the geodesic path between the initial and final equi- 
librium states in advance and then apply the stepwise algorithm along the 
geodesic path. It is worth noting that in simulated annealing |J the above 
algorithm can not be directly applied since the annealed system's temper- 
ature and specific heat are known but statistically. We nevertheless think 
that our results could be adapted to computer simulated statistical systems. 

A few words are needed to interprete the condition of fixed thermody- 
namic time when looking for the optimum control. It may often happen that 
it is not the clock-time of a given thermodynamic process that is the eco- 
nomically or technologically relevant quantity. Rather than clock-time, the 
number of iterated technological steps might better characterize the amount 
of relevant facilities (computer capacity, special materials, or just money) 
that can be consumed to bring the system from its initial state into a pre- 
scribed final one. In such cases the clock-time period of the process is of less 
interest to be fixed in advance. The step number is rather to be fixed. The 
condition of fixed-in-advance step number is equivalent to the condition of 
fixed-in-advance thermodynamic time, as shown by Eq. (f|0]). It may eventu- 
ally happen that the thermodynamic time (proportional to the step number) 
takes the place of the clock-time. The overall clock-time needed to perform 
a single cooling step might not be dominated by the clock-time of the relax- 
ation but by the clock-time of the technological adjustment before and after 
the relaxation. In this case, the gross time of the stepwise process will be 
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proportional to its thermodynamic time. Actually, the fixed-in-advance ther- 
modynamic time becomes the relevant condition in designing the finite-time 
thermodynamic process and, consequently, minimum entropy production will 
be achieved along the geodesic path in thermodynamic state space. 

6 Conclusion 

We have extended the notion of the relaxation time r from one dimension to 
an arbitrary number of dimensions. Thus the thermodynamic time can be 
deduced for multi-dimensional thermodynamic processes, too, in a straight- 
forward way. As a consequence, the optimum control of a thermodynamic 
process based on the thermodynamic time scale can be achieved with any 
number of variables to be controlled. We have seen that the optimum path 
should be followed as close as possible by the process in hand as the geodesic 
is uniquely marked out between the initial and the desired final states. In 
addition, the proposed stepwise driving force for the process can be arranged 
by measuring ordinary thermodynamic variables e.g. temperature, concen- 
tration etc. 

A remarkable result of the present analysis is that the path of mini- 
mum irreversible entropy production in fixed thermodynamic time becomes 
a geodesic of the Riemannian space of thermodynamic states introduced a 
decade ago. Thus far Riemannian structure was seen only in its influence 
in infinitesimal neighbourhoods, e.g., on thermodynamic fluctuations. Now 
we have shown within the context of standard irreversible thermodynamics 
that some finite time thermodynamic processes, when optimum, follow the 
shortest paths of the Riemannian space. The result is valid in an arbitrary 
number of dimensions and proves for the first time the distinguished role of 
geodesies in driven irreversible processes. Our result is based on the usual 
expansion of the irreversible entropy production up to the second order in 
speeds. Of course, the given interpretation for geodesies is restricted to not 
too high speeds and rates. (This, however, would not at all restrict the ul- 
timate value of the interpretation. Recall the analogous case in space-time 
geometry: it is not the too large objects but the small enough ones that will 
follow geodesies.) These results may turn out to form a basis for amplifying 
the benefits of finite-time thermodynamics both in theory and in practice. 

After the completion of our work we became aware of a related report by 
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Spirkl and Ries [ j"0"| . Whereas they do not use Riemann-geometrical terms 
there seem to be significant coincidences between our results. 
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